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f*^; Herein we show that the historical records of mid-latitude auroras from 1700 to 1966 present oscillations with periods of about 9, 
CN| 10-11, 20-21, 30 and 60 years. The same frequencies are found in proxy and instrumental global surface temperature records since 
^^1650 and 1850, respectively and in several planetary and solar records. Thus, the aurora records reveal a physical link between 
C3 climate change and astronomical oscillations. Likely, there exists a modulation of the cosmic ray flux reaching the Earth and/or 
\! of the electric properties of the ionosphere. The latter, in turn, have the potentiality of modulating the global cloud cover that 
ultimately drives the climate oscillations through albedo oscillations. In particular, a quasi 60-year large cycle is quite evident 
since 1650 in all climate and astronomical records herein studied, which also include an historical record of meteorite fall in China 
from 619 to 1943. These findings support the thesis that climate oscillations have an astronomical origin. We show that a harmonic 
,—i constituent model based on the major astronomical frequencies revealed in the aurora records is able to forecast with a reasonable 
accuracy the decadal and multidecadal temperature oscillations from 1950 to 2010 using the temperature data before 1950, and 
vice versa. The existence of a natural 60-year modulation of the global surface temperature induced by astronomical mechanisms, 
r~ | by alone, would imply that at least 60-70% of the warming observed since 1970 has been naturally induced. Moreover, the climate 

Q_imay stay approximately stable during the next decades because the 60-year cycle has entered in its cooling phase. 
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1. Introduction 

Since ancient times people have claimed that climate and 
weather changes are related to cyclical astronomical phenom- 
ena linked to the orbits of the Sun, the Moon and the planets 
( Ptolemv, 2nd centum iMa'sarl 1886c iKeplen, Il606t ISwerdlowl 
1998; Iyengar, 2009). Because of this conviction the ancient 
astronome r s deve loped calendars that contain several cycles 
jlAslaksenl 1 1999b as well as the well-known annual cycle. 
During the last 70 years, numerous scientific evidences appear 
to have cor r oborat ed that ancient conviction. For example, 
MilankovicJ dl94lb theorized that variations in eccentricity, 
axial tilt and precession of the orbit of the Earth determine 
climate patterns such as the 100,000 year ice age cycles of 
the Quaternary glaciation. Milankovitch's theory fits the data 
very well, over the past million years, in particular if the 
temporal rate of change of global ice volume is co nsidered 
( Roe et all 2006b. More recent l y, a number of aut hors dShavivl 
2003? Shaviv and Veizer . 20031: Svensmarki 2007 ) have shown 
that the cosmic -ray flux records well correlate with the warm 
and ice periods of the Phanerozoic during the last 600 mil- 
lion years: in this case the cosmic-ray flux oscillations are 
believed to be due to the changing galactic environment of 



the solar system, as it crosses the spiral arms of the Milky 
Way. Over millennial and secular time scales several authors 
have found that changes in sunspot number and cosmogenic 
isotope produ c tions well corre la te wi t h climate changes 
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2008 F iRaspopov et all |2008[ lEichler et all 120091 ISoonl 12009 
iMeehl et aiT l2009; Sca fettal l2009b . Moreover, instrumental 
global surface records since 1850 appear to be characterized 
by a set of frequencies that can be associated to the Moon 
(9.1-year period) and to the motion of the Sun relative to the 
barycent er of the solar sy stem (about 10.5, 20, 30 and 60 year 
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In this paper, we study the historical mid-latitude aurora 
records since 1700 (IKfivskv and Peimllll988tlSilvermanill992l) 
and show that these records share the same set of frequencies 
that characterize the climate system as well as the natural oscil- 
lations of the solar system. This finding reveals the existence 
of a clearer phy sical mechanism, missing in our previous study 
dScafettal. 12010b). that could link the astronomical cycles to cli- 
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mate oscillations. The major implication of this paper is that 
there exists an astronomical harmonic modulation of the elec- 



cover and, therefore, the terrestrial albedo ( 


Svensmark 1998 


2007: Svensmark et al.. 20091 Kirkbv. 2007: 


Tinslevl 2008b. 



This research would also support the development of a novel 
astronomically-based theory of climate change that may cred- 
ibly compete and likely substitute the current mainstream an- 
thropogenic g lobal warming theory (AGWT) advocated by the 



IPCd ( 120071) . In fact, the AGWT advocates claim that as- 
tronomical forcings of the climate are almost negligible and 
that the climate variations are induced by some still poorly 
understood and modeled internal chaotic dynamics of the cli- 
mate system, and by trends in greenhouse gases (GHG) (mostly 
CO2 and CH4) and aerosol re cords. More precise ly, global 
surface temperature has risen dBrohan et all 2006 1 by about 
0.8 °C and 0.5 °C since 1900 an d 1970, respec t ively The 
IPCd d2007l) and other researchers (Lean and Rindi E0O8 ) have 



claimed that more than 90% of the observed warming since 
1900 and practically 100% of the observed warming since 
1970 have had an anthropogenic cause. While Lean and 
Rind's methodology that the climate responds linearly with 
the forcings can be easily questio ned by noting th at the heat 
capacity of the Earth is not zero (IScafettai 120091) . the IPCC 
claims derive from figures 9.5 and 9.6 in the IPCC report 
(AR4-WG1, 2007) showing, by means of professional cli- 
mate general circulation models (GCMs), that natural forc- 
ings alone (volcano and solar irradiance) would have caused 
a cooling since 1970: thus, the observed post- 1970 warm- 
ing has been interpreted as being induced by human ac- 
tivity alone. However, the very large uncertainty in the 
aerosol forci ngs and of the clima te sen sitivity to GHG changes 
dlPCd 120071 iKnutti and Hegerl 120081: [Rockstrom etail 120091: 
iLindzen and Choii 1201 It ISpencer and BraswellL 1201 II) . and 
the current very poor modeling of the wa ter vapor feedback 
( Solomon et al. , 20 1 Oh . of the cloud system ( Lauer et al. , 2010). 
of the ocean dynamics and of the biosphere question the robust- 
ness of the current GCM s for properly interpret ing and recon- 
structing the real climatedds o and Singerl l2009l) . 

On the contrary, if the climate system is mostly character- 
ized by a specific set of harmonics, it may be possible to par- 
tially reconstruct and forecast it in the same way in which ocean 
tides are currently predicted, that is, by means of harmonic con 
stituent models based on astron omical cycles ( Thomson! 1881 



Fischer, 1912 : Doodsonl 1921). A harmonic constituent model 



is just a superposition of several harmonic terms of the type 



N 



F(t) = A + ^ Ai cos(<y,-f + 



(1) 



The frequencies are deduced from astronomical theories. 
The amplitude A, and the phase of each harmonic are em- 
pirically determined using regression on an adequate sample of 
observations, and then the model is used to forecast future sce- 
narios. Currently, in most US costal loc ations tidal f orecast is 
made with 35-40 harmonic constituents (Ehret, 2008). 



In the following, we show that typical energy balance models 
and general circulation models can be mathematically reduced 
to harmonic models in first approximation, and propose that 
the climate oscillations too can be approximately reconstructed 
and forecasted by using a planetary harmonic constituent model 
philosophically equivalent to Eq. [TJbased on astronomical cy- 
cles. 



2. A possible link between mid-latitude auroras and the 
cloud system 



In this paper, we postulate that the annual frequency occur- 
rence of mid-latitude aurora events is a measure of the level of 
electrification of the global ionosphere, which is most l y reg- 



ulated by incom ing cosmic ray flux variations (iKirkbvi 12007 



Svensmarld 120071) . When the ionosphere is highly ionized by 
cosmic rays, large auroras would more likely form at the mid- 
latitudes. This phenomenon would occur because when the up- 
per atmosphere is highly ionized, it would also be electrically 
quite sensitive to large solar wind particle fluxes and favor the 
formation of extended mid-latitude auroras. In fact, higher ion- 
ization of the atmosphere would mostly occur when the magne- 
tosphere is weaker and cosmic ray as well as solar wind parti- 
cles, can more easily reach the mid-latitudes. Then, the level of 
atmospheric ionization and of the global electri c circuit of the 



atmosphere should regulate the cloud system dKirkbv , 2007 



Svensmar k, 2007; Tinslev, 2008). If the above theory is cor- 
rect, the frequencies of the mid-latitude aurora records should 
be present in the climate records too. 

Indeed, cloud-related climatic effects can largely d ominate 
other mechanisms such a s CO2 and CH4 GHG forcing (Kirkbv, 
2007; Svensm ark! l2007h . For example, in the past billion years 
the Earth experienced severe glaciations despite the fact that 
the CO2 concentration was at least 10 times higher than to- 
day ( 4000-6000 ppmv against the actual 390 ppmv) dHavdenl 
2007b . In add ition to major continental d rifts and other geo- 
logical events (ICourtillot and Rennd. 12003). it has also been ar- 
gued that the high cosmic ray incoming flux that occurred dur- 
ing known major glaciations could have i ncreased the cloudi- 
ness of the Earth causing a global cooling dShavivi l2003i [2008 ; 
Shaviv and Veizeii 120031; ISvensmarkll2007l) . 

The cloud system controls a large part of the terrestrial 
albedo and regulates the amount of total solar irradiance reach- 
ing the Earth's surface. The solar irradiance reaching the 
Earth's surface warms the ocean and the land. A small astro- 
nomical modulation of the terrestrial albedo through the cloud 
syste m can greatly i ncrease the clim ate sensitivity to solar forc- 



ing dShavivi 120081 IScafettai l2009h . Evidently, if the current 
GCMs are missing an important forcing of the cloud system, 
they would poorly reconstruct and, potentially, severely mis- 
interpret climate variations at multiple temporal scales. In- 
deed, a correlation between galactic cosmic ray fluxes and 
high-altitude, mid-altitude and low clouds have been found 
dSvensmarkl 120071 iRohset allbOlOtrLaken et al.Ll2oTob . 
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Figure 1: Global surface temperature anomaly (Brohan et al.| 2006). The fig- 
ure also shows the quadratic fit upward trend of the temperature (black). 



3. A 60-year cycle in mid-latitude aurora and climate 
records 

Figure 1 sh o ws th e global surface temperature (HadCRUT3) 
dBrohan et all 120061) from 1850 to 2010 (monthly sampled). 
Global surface temperature records (land, ocean, north hemi- 
sphere, south hemisphere) have been found to be char acterized 
by a cl ear and large quasi 60-year cyclical modulation (Scafetta, 
l2010bl) . In fact, the following 30-year trends are evident in 
the record: 1850-1880, warming; 1880-1910, cooling; 1910- 
1940, warming; 1940-1970, cooling; 1970-2000, warming; and 
a small cooling since 2000 that may last until 2030-2040. 

Figure 2A shows the global surface temperature record de- 
trended of its upward trend and smoothed with a 8 -year moving 
average that highlights its 60-year cyclical modulation with a 
peak-to-trough amplitude of about 0.3-0.4 °C. 

Figure 2B shows the annual frequencies of mid-latitude au- 
roras obtained from the supplement of the catalogue of mid- 
latitude auroras <55N from 1700 to 1900. This record contains 
the historical aurora obs ervations reported mostl y in central Eu- 
rope since 1000 AD dKfivskv and Peimll, 1 19881) . Before 1700, 
the record is largely incomplete and the data are not depicted in 
the figure. Figure 2B also depicts the catalog of the aurora ob- 
servations in the Faroes Islands from 1872 to 1966. Despite the 
fact that the Faro es' record refers to a northern region (62N), 



Silvermanl (1 19921) noted that it appears to have physical prop- 
erties compatible with the Krivsky and Pejml's record. So, it 
may be reasonable to combine the two catalogs covering 266 
years from 1700 to 1966. The present author does not know 
why these two records stop in 1900 and 1966 respectively. Au- 
rora events surely have be collected but not with the same naked 
eye methodologies used in the past. The most recent electron- 
ically based records could show technical heterogeneities with 
the historical ones such that the old and recent records cannot be 
simply combined. The aurora records are flipped up-down and 
superimposed to the same 60-year periodic sinusoidal function 
depicting the 60-year temperature cycles in Figure 2A. Mid- 
latitude auroras are more likely to occur during cooler multi- 
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Figure 2: [A] The 60-year cyclical modulation of the global surface temper- 
ature obtained by detrending this record of its upward trend shown in Figure 
1. The temperature record has been filtered with a 8-year moving average. 
Note that detrending a linear or parabolic trend does not significantly deform 
a 60-year wave on a 160-year record, which contains about 2.5 of these cy- 
cles, because first and second order polynomials are sufficiently orthogonal to a 
record containing at least two full cycles. On the contrary, detrending higher or- 
der polynomials would deform a 60-year modulation on a 160-year record and 
would be inappropriate. [B] Aurora records from the "Cat alogue of Polar Au- 
rora < 55N in the Period 1000-1900" from 1700 to 1900 l lKrivskv and Peimll 
1988). Figure 2B also depicts the catalog referring to the aurora observa- 
tions from the Faroes Islands from 1872 to 1966. Both temperature and aurora 
records show a synchronized 60-year cyclical modulation as proven by the fact 
that the 60-year periodic harmonic function superimposed to both records is the 
same. 



decadal periods. 

An approximate 60-year modulation in the au- 
rora record from 1000 to 1900 has been observed by 

once that the original 



Charvatova-Jakubcova et al.l (119881) . 
record before 1700 is adjusted with wha t is th ere called a 
"civilization" factor coefficient. iKomitovl ( 120091) also found 
a correspondent quasi 60-year cycle in both aurora and 
Greenland and Antarctic U) Be concentration records during the 
period 1700-1900. Periodic patterns and long-term variations 
in the aur ora records have b een observed for centuries since the 
work by de Mairan) dl733 ). De Mairan found s everal reprises 
and returns in the aurora activity. Siscoe (1980) summarized a 
number of the studies that searched for periodicities in aurora 
occurrence. For exam p le, in 1831 Hansteen inferred a period 



of 95 years; lOlmsted (1856) found distinctive evidences for 
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Figure 3: [A] Twenty-year moving average of the tree-ring chronologies from 
Pinus Flexilis in California and Albertain: th is record is used as a proxy f or re- 
constructing the Pacific Decadal Oscillation (MacDonald and Case, 2005). [B] 
Record of G. Bulloides abundance variat ions (1-mm intervals) from 1650 to 
1990 A.D. (black line) iBlack et all l 19991; this is a proxy for the Atlantic vari- 
ability since 1650. [C] Five-yea r running average of the Ind ian summer mon- 
soon rainfall from 1813 to 1998 lAgnihotri and Dutta, 2003). All three records 
show clear 60-year cyclical modulations that are (positively or negatively) well 
correlated to the 60-year cycles of the global surface temperature and the aurora 
records. The records are best fit with sinusoidal functions that give a statistical 
error about the 60-year period of ±4 years. 



60 to 65 years as the intervals between great auroras; and 
Wolf and Fritz preferred a period of about 55 years. Another 
work also indicated periods in the region of 80-90 years 
dFevnman and Fougerel 1 1 9841) . and even large periods as long 
as 200 or 400 years. Large cycles with periods of about 60, 
80-90 years and longer bi- and multi-secular cycles are the 
most commonly reported. These frequencies correspond to the 
well-known Gleissberg and Suess solar cycle frequency band s 
dSuesslI 1 9801: iKfivskv and Peimllll988l;IOgurtsov et alfl2002h . 

As already discussed in LScafettal d2010 ab), quasi 60-year 
cycles are found in s everal climatic from several region s of the 
Earth (fKlvashtorinl l200lt iKl vashtorin and Lvubushin, 2007; 
Klvashtorin et all 120091 : iLe Mouel et all 120081; ICamuffo et ail 
2010 1: lAgnihotri et all 120021: lAgnihotri and Duttal 12003 



Sinha et all 120051: iGoswami et al.lT2006l: lYadava and Ramesh] 
2007!: iMazzarella and Scafettal 1201 ll Ijevreieva et all 12 008 ) 
and astronomical/solar records (lYu et allll983HPatterson et al 



clearly suggest the existence in the climate of a natural 60-year 
cycle synchronized to a correspondent solar/astronomical cy- 
cle. For example, Figure 3 depicts three multi-secular climate 
records from three independent regions of the globe showing 
multiple quasi-60 year large oscillations since 1650. Figure 3 A 
depicts a record obtained from the tree-ring chronologies from 
Pinus Flexilis in California and Albertain: the best sinusoidal 
fit gives a period of T = 61.5+4 years. This record is used 
as a proxy for reconst r uctin g the Pacific Decadal Oscillation 
dMacDonald and CaseL 120051) . Figure 3B depicts the G. 
Bulloides abundance variation record found in the Cariaco 
Basis sediments in the Caribbean sea since 1650 dBlack et all 
1999b : the best sinusoidal regression fit gives a period of 
T = 59.5 + 4 years: a quasi-60 year cycle has been found in 
the Caribbean sea fo r millennia during the entire Holocene 
( Knudsen et al. . 19881) . This record is an indicator of the trade 



2004; Qg urtsov et a l.. 2002; Robe rts et all l2007). These results 



wind strength in the tropical Atlantic Ocean and of the North 
Atlantic Ocean atmosphere variability. Periods of high G. 
Bulloides abundance correlate well with periods of reduced 
solar output (the well-known Oort, Wolf, Sporer, Maunder, 
Dalton mi nima): a fact that s uggests that these cycles are solar 
induced (Blac k et all Il999h . Figure 3C depicts the Indian 
summer monsoon rainfall record from 1813 to 1998 years 
that also shows prominent quasi 60-year cycles. The Indian 
summer monsoon rainfall record, together with east equatorial 
and Chinese monsoons, clea rly manifest a solar variab ility 
signature for several centuries (lAgnihotri and Dutta 2003): the 
best sinusoidal fit of the Indian monsoon rainfall record gives 
a period of T = 62 + 4 years. The cycles depicted in Figures 2 
and 3 are well synchronized to the quasi 60-year modulation of 
the global temperature observed since 1850. 

Consequently, the very good correspondence between the 
two 60-year periods 1880-1940 and 1940-2000, which are ob- 
served in all (North a nd South , Ocean and Land) global surface 
temperature records (Scafettal 2010allb ). is unlikely just a co- 
incidental red-noise pattern. With a high confidence level, a 
natural 60-year periodic modulation exists in the climate sys- 
tem because these cycles are observed in a large number of 
natural multi-secular records although this cycle may not ap- 
pear always evident because of errors in the data and because 
of other superimposed patterns of different physical origin (for 
example, volcano effects) and other natural cycles that could 
induce interference patterns. 

4. Spectral analysis of aurora and climate data, and joint 
power statistical tests 

To determine whether faster cycles are equally common in 
the aurora and temperature records we analyze the global sur- 
face temperature and the auroras re cords by using the maximum 
entropy method dGhil et all [2002) because with a proper high 
number of poles this technique resolves the very low frequency 
band of the spectru m much better than Fourier or periodogram 
based techniques dPriestlv et al. . 2001). Note that in our case 
MEM needs to be used with a very large number of poles (up 
to half of the length N of the record) because we are inter- 
ested in resolving the very low frequency range of the spectrum 
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Figure 4: [A] Power spectrum of the global surface temperature [A], which 
covers the period 1850-2010. [B,C,D] power spectra of three auroras records, 
which cover the periods 1700-1966, 1700-1880 and 1872-1966 respectively. 
The numbers #1, #2, #3, #4 and #5 indicate the periods at about 9, 10-11, 20, 
30 and 60 years. These peaks are statistically significant against red-noise back- 
ground tests at different confidence levels (median, % 90, %95, %99), which 
are indicated with the dot lines. 



where the dacadal and multidecadal periodicities are located, as 
simple computer experiments with synthetic harmonic records 
would easily prove. A number of p oles between N/5 and N/2 



(at most) is correct in most cases (Ul rvch and BishopL 11975 



Couilillot et al., 1977) 



Figure 4A shows the power spectrum of the temperature that 
covers the period 1850-2009 . This curve corresponds to those 
analyzed in IScafettal (bOlObl) . where additional alternative sta- 
tistical tests to check the robustness of the results are conducted. 
Four major decadal and multidecadal cycles are seen at about 
9, 10-11, 20-21 and 60-year periods. These cycles are indicated 
with the number #1, #2, #3 and #5. With the exception of the 
smaller peak #4 (about 30 years), which may not be signifi- 
cant, all other four frequencies have a 99% confidence against 
red noise background. Faster cycles or less important cycles 
are not studied here. Traditional periodogram analysis (not re- 
ported in this figure) produces the same major four peaks de- 
tected by MEM, but with a larger uncertainty, and confirms that 
the global surface temperature is characterized by at least the 
above four major frequency peaks. 

Figure 4B shows the power spectrum of the superposition of 
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Figure 5: Joint power statistics l lSturrock et aD.l2005l) between the power spec- 
tra of the 1850-2010 global temperature (Figure 4A) and the 1700-1900 aurora 
record (Figure 4C). The joint peaks at periods of about 9, 10.5, 20 and 60 years 
are evident and have a 99% confidence level against red-noise background. 



the catalogues of mid-latitude auroras merged with the Faroes' 
auroras: this combined record covers the period 1700-1966. 
Again, five peaks cycles are clearly seen at about 9, 10-11, 
20-21, 30 and 60 years. Figure 4C shows the power spectrum 
of a portion of the catalogue of mid-latitude aurora covering 
the period 1700-1880, which presents a minimum overlap with 
the temperature data and with the Faroes' catalog. Again, five 
peaks are seen at about 9, 10-11, 20-21, 30 and 60 years. Figure 
4D shows the power spectrum of the Faroes' catalog that cov- 
ers the period 1872-1966. In this case, only 3 cycles are clearly 
seen: cycle #2 (10-11 years); cycle #3 (20-21 years); and cycle 
#5 (about 60 years). Note that cycle #1 (about 9 years) is not 
as strong as the one found in the temperature record: we will 
discuss this finding later. 

Table 1 reports the measured five periods and their average 
values. The frequencies of the temperature and of the AAR 
columns are compatible within their error of measure. For the 
5 period couples reported in Table 1, the reduced x 2 is xl = 
0.18 with 5 degrees of freedom and the coherence probability 
isP 5 (x 2 >xl) >96%. 

Figure 5 further supports the result that the aurora records 
and the global temperature record share a similar set of fre- 
quency. In this case, we use the joint power statistic function 



(Sturrock et al., 2005) which confirms that the aurora record 



and the global temperature record share a common set of major 
cycles at about 9, 10.5, 20, and 60-year periods (cycles #1, #2, 
#3 and #5) with a 99% confidence. If two power spectra share a 
compatible frequency peak, also their joint power statistic func- 
tion, which is approximately given by ^jS \(f) * S2(f), would 
show a significant peak at the same frequency. Other common 
cycles may be present as well, but they are less evident. 



5. A possible astronomical planetary origin of the aurora 
and temperature oscillations 

The observed aurora and climate cycles suggests that the 
solar system planetary orbital dynamics m ay be their first 
physi c al cause. I ndeed , some studie s dWolil 1 18591: ISchusteij, 



1911 



1987 



Bendandi, 1931; Jose, 119651: 'airbridge and Shirle 
Landscheidl 1 19881 1 1 999 t ICharvatoval 1 1 990i l2000l 1200 



Charvatova and Stfestfld 120041: iGrandpierrel Il996t iMackevl 
2007tlHungll2007tlWilson et alJj2008t|Sidorenkov and Wilson! 



2009t IPerrvman and Schulze-Hartungl uOlOt) have suggested 



that solar variation can be partially driven by the planets 
through gravitational spin-orbit coupling mechanisms and grav- 
itational tides. 

A full theory that would physically explain how the solar 
wobbling or the planetary tides may influence solar activity 
has not been fully developed yet. However, preliminary stud- 
ies suggest that planetary gravity may incre ase nuclear rate 
( Grandpierrel 1 19961 I Wolff and Patronet l2010h by favoring the 
movement of fresh fuel into the solar core. The proposed mech- 
anisms would likely produce the major frequencies herein dis- 
cussed because it is based on the study of the wobbling of Su n 
around the solar system bary center as done in Scafetta ( 2010bl) . 

Moreover, solar wobbling could be physically relevant rela- 
tive to the incoming cosmic ray flux. In fact, because the inner 
part of the solar system is almost gravitationally locked to the 
Sun, and the Earth's orbit too wobbles around the solar system's 
barycenter almost in the same way as the Sun does. Indeed, an 
observer from the Earth does not see the solar wobbling while 
an observer external to the solar system would see the Earth's 
orbit wobbling almost in the same way as the Sun does. Thus, 
relative to the Sun or to the Earth's orbit, the incoming cosmic 
ray flux is wobbling wit h the same freq uencies found in the so- 
lar barycentric motion (IScafet ta, 2010b). 

The major planetary resonances of the solar system related 
to Jupiter and Saturn are the following: Jupiter's sidereal pe- 
riod 11.862 years; Saturn's sidereal period 29.457 years; the 
opposition-synodic period of Jupiter and Saturn (9.93 year), 
which oscillates between 9.5 and 10.5 years due to the eccen- 
tricity of their orbits; about 19.84 years, the synodic period of 
Jupiter and Saturn; about 59.6 years, the repetition of the com- 
bined orbits of Jupiter and Saturn; and a spring tidal beat fre- 
quency (~61 year). Moreover, there is the 1 1-year Schwabe so- 
lar cycle and the 22-year Hale magnetic solar cycle which per- 
haps are indirectly associated to the two major Jupiter/Saturn 
tides because the 11 -year cycle appears constrained between 
the 9 .3-year spring Jupiter/Saturn tide and the 1 1 .8-year Jupiter 
tide IWilsorJ ( 1987 ). Thus, three major solar and solar system 
oscillations should be expected at periods of about 10-11, 20- 
22 and 59-63 years. Many other resonance frequ encies related 



to all planets have been found in the Sun s motion (Buc ha et al 



119851: ICharvatova-Jakubcova et all 11988), but they are not dis- 
cussed here. 

An additional large 9.1 -year temperature cycle doe s not ap- 
pear among the planetary resonances and IScafettal d2010b) 
found that it may be a solar/lunar tidal cycle. In fact, a 9.1 -year 
periodicity is between the period of the recession of the line of 



lunar apsides, about 8.85 years, and half of the period of pre- 
cession of the luni-solar nodes, about 9.3 years (the luni-solar 
nodal period is 18.6 years). The 9. 1-year cycle is also about half 
of the 18.03-year Saros eclipse cycle. Note that every almost 
9 tropical years the Sun, the Earth and the Moon are aligned 
at the same angle relative to the Earth's surface. Because two 
opposite tides are simultaneously present, the two alignments 
Sun-Moon-Earth and Sun-Earth-Moon are physically equiva- 
lent. Consequently, a quasi 18-year astronomical solar/lunar 
cycle is associated to a quasi 9-year tidal cycle dynamics. 

The results depicted in Figures 4B, 4C and 4D sug- 
gest that the aurora events are synchronized to the above 
astronomical periodicities. About the two decadal cy- 
cle in the climate records, several studies have attempted 
to interpret it either as a lunar influence on climate 
via the 18.6-year luni-solar nodal or the nutation cycle, 
or as an infl u ence o f the 22-year Hal e solar magnetic 



cycle dCurriel 1 19841 iKeeling and Whorfl I2OO0I ICamuffc 



200 it iMcK innell and Crawford! l2007t iMunk and Bills] . 12007 
Hovt and Schattent 1997b . However, on a global scale the most 



prominent bi-decadal cycle has a period of about 20-21 years 
as shown above in Figures 4 and 5, while a distinct 18.6 year 
lunar cycle or a 22-year solar cycle are not clearly visible in the 
adopted global temperature records. Perhaps, a 18.6-year nu- 
tation cycle may become particularly significant only in some 
regions, such as the polar ones. Because the global data herein 
analyzed appear to show a quasi 20-21 year cycle more than 
a sharp 18.6 or a 22-year cycle, we claim that a quasi 20-21 
year cycle is w hat may characte rizes the well known bi-decadal 
climatic cycle (Scafetta, l2010bh . 



6. Further analysis of Jupiter and Saturn's orbital dynam- 
ics and of their induced solar system cycles 

Herein we further analyze how the orbital eccentricities and 
the periods of Jupiter and Saturn can produce major 10, 20 
and 60-year oscillations in the solar system. Figure 6 shows 
the dates of all conjunctions of Jupiter and Saturn occurred 
from 1650 to 2010 relative to the Sun. The orbital data are 
obtained using the NASA Jet Propulsion Laboratory Develop- 
mental Ephemeris. The angular separation between Jupiter and 
Saturn relative to the Sun at the conjunction day is of the or- 
der of 1", or below. Table 2 reports the dates, the heliocentric 
coordinates and the distance from the Sun of the conjunctions 
of Jupiter and Saturn. The coordinates are chosen to be the 
right ascension and declination, which refer to the equatorial 
plane of the Earth. Thus, the declination measures the direction 
of the gravitational and magnetic forces of Jupiter and Saturn 
relative to the daily average orientation of the terrestrial mag- 
netosphere. The figure shows that the conjunctions since 1650 
can be separated into three groups. A conjunction in the Group 
I is followed by a conjunction in the Group II at an angle of 
about 360° * 19.86/29.46 = 242.7° after about 20 years, which 
is followed by a conjunction in the Group III again at about 
242.7°, after about 20 years. Every 3rd alignment (about 60 
years) is about 8.1 degrees of the original starting point: the 
configuration repeats about every 900 years. A similar diagram 
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Km and 1600 Km respectively. The results about the tidal elon- 
gation indicate that Group I alignments occur when Jupiter and 
Saturn occupy the closest positions to the Sun and their com- 
bined tidal effect is the largest. This returning pattern gives 
origin to a 60-year physical tidal cycle on the Sun and in the 
heliosphere in proximity of the Earth's orbit. A similar 60- year 
periodic pattern would be generated by the magnetic fields of 
Jupiter and Saturn that influence the heliosphere too. 

To show the varying effect of the tidal combination of the J/S 
synodic cycle and Jupiter eccentricity cycle we observe that: 
a) the tidal oscillation produced by Jupiter and Saturn orbits at 
1 AU have an amplitude of about 224 km and 12.3 Km, re- 
spectively, which are calculated as half of the difference be- 
tween the tidal elongation at the perihelion and at the aphe- 
lion; b) J/S synodic cycle produces an average tide amplitude of 
about 38.4 Km, which is the half amplitude of the average Sat- 
urn induced tidal elongation. Jupiter's perihelion occurred on 
20/May/1999 (n = 1999.38), Saturn's perihelion occurred on 
27/Jul/2003 (r 2 = 2003.57) and the J/S conjunction occurred 
on 23/Jun/2000 (r 3 = 2000.48). Thus, the average J/S tidal 
anomaly effect approximately oscillates as 



Figure 6: Dates and positions of the conjunctions between Jupiter and Saturn 
from 1650 to 2010, as reported in Table 2. The figure depicts the orbits of the 
Earth (E), of Jupiter (J) and of Saturn (S). 'P' is the position of the perihelion of 
the Earth. The figure emphasizes the tri-synodic cycle of Jupiter and Saturn. It 
takes about 60 years for Jupiter and Saturn to reach the same relative alignment 
around the sun. Two J/S conjunctions occurred in 1146 and 1206 are added. 
The latter two J/S conjunctions reveal a 800-850 year cycle in the J/S conjunc- 
tion circulation as seen from the Earth and noted by Kepler who was using the 
tropical orbital period of the planets. 



was prepared by Kepler ( 1606) to interpret the Star of Betlehem 
as a J/S conjunction ( Kemp, 2009), but Kepler used the tropical 
orbital periods that yield an angular separation of about 3 de- 
gree: and the configuration would repeat approximately every 
800 years, as depicted in Figure 6. 

The last column of Table 2 reports the approximate magni- 
tude in Km of the tidal elongation produced by the combined 
effect of Jupiter and Saturn at 1 AU from the Sun, which is the 
Sun-Earth average distance. The tabulated spring tidal elonga- 
tion values at 1 AU from the Su n are calculat ed by using the 
following approximated formula (Taylor 2005 ): 



3 M, 



a+ - a- = — 



uE 



3 M s 
+ - 



^SuE 



(2) 



where: a+ and a_ are, respectively, the maximum and minimum 
radius of the equipotential surface spheroid produced by the 
gravitational field of the Sun plus the tidal potentials of Jupiter 
and Saturn; M Su = 333000, Mj = 317.9 and M s = 95.18 are 
the masses of the Sun, Jupiter and Saturn in Earth's masses; 
Rs,ij and Rs„s are the distances of the Sun from Jupiter and 
Saturn, respectively; and Rs u e = 1 AU is the average Sun-Earth 
distance. 

Table 2 shows that the tidal elongations at the J/S conjunc- 
tions occurring in Group I are significantly larger (about 1800 
Km) than those occurring in the other two groups: about 1400 



D(f) = 224 cos 



2n (f-Ti) 



11.86 



(3) 



38.4+ 12.3 cos 



2tt (t - T 2 ) 



29.46 



cos 



2k (t - T 3 ) 



9.93 



where it is assumed the superposition of two waves given by 
the tidal sidereal orbit of Jupiter and the spring tidal oscillation 
of Jupiter and Saturn. 

Figure 7 A depicts the function D(t) for the period 1650-2050. 
The quasi 60-year beat modulation is evident, and it is in good 
phase with the 60-year modulations observed in Figures 2 and 
3. Note, for example, that the function D(t) presents minimum 
amplitudes during periods of temperature minima (1850-1860, 
1900-1920, 1960-1980) and maximum amplitudes during pe- 
riods of temperature maxima (1870-1890, 1930-1950, 1990- 
2010). 

Figure 7B depicts the speed of Sun relative to the barycenter 
of the solar system where the 20-year and 60-year modulation 
of the solar orbit is also evident. Note that the 60-year mod- 
ulation of the solar speed corresponds exactly to the 60-year 
modulation of the auroras and of the climate depicted in Figure 
2. 

The existence of an astronomical 60-year cycle that is in- 
fluencing the inner part of the solar system is furthe r proven 
by m eteorite fall records, as extensively studied in lYu et al 
( 1983.1). It is well known that Jupiter and Satur n can regulate 



the collision rate of small bodies with the Earth (Bennett et al 



20101) . Figure 8 shows our analysis of a number of historically 
recorded meteorites falls in China from AD 619 to 1943. Al- 
though this record may be incomplete and should be corrected 
by some social parameters, it shows cycles at about 10.5, 30 
and 60 year periods. This result would stress the importance 
of the J/S 60-year cycle in driving a major astronomical oscil- 
lation of the solar system. The power spectrum in the insert of 
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Figure 8: N umber of histori cally recorded meteorite fall in China from AD 
619 to 1943 Yu et al. 1 1983). Note that the gradual multi-secular increase of 
witnessed meteorite falls is likely due to social parameters effects. In fact, in 
China, the typography (book-printing) has been known since AD. 1041 while 
only hand written records are available before, so the records are very scarce 
in the first millennium in comparison with the greater number of records in 
the second half of the second millennium. However, this trending should not 
significantly alter the faster decadal and multidecadal modulation, which should 
be physical meaningful. The insert shows the power spectrum of the record with 
the Maximum Entropy Method (50 poles) and the Multi Taper Method. Cycles 
at 10.5, 30 and 60 years are clearly visible. 



Figure 7: [A] The Jupiter-Saturn tidal anomaly as deduced from Eq.[3]at 1 AU 
from the Sun. Note the 60-year beat modulation. This 60-year modulation is in 
good phase with the 61 -year cycle of the temperature and of the aurora records 
since 1700 depicted in Figures 2 and 3: the largest magnitudes correspond to 
peaks of the 61-year cycles. [B] Figure B depicts the speed of Sun relative to the 
solar system barycenter (SSSB) (Scafetta, 2010b). The 20-year (solid) and 60- 
year (dash) oscillation is evident. The dash curve is a moving average filtering. 
Note that the 60-year modulation observed in both figures corresponds exactly 
to the 60-year modulation of the auroras and of the climate depicted in Figure 
2. 



Figure 8 also suggests the astronomical origin of a significant 
~10.5-year cycle, which is also found in the temperature record 
(see Table 1) and in t he alignments of the tide-raising plan 
ets within ten degrees and Anderson! 1 1 975 : Buc ha et al 
ll985HGrandpierrel Il996b . 



Quasi-millennial solar c ycles are evident in several cosmo 
genie isotope production s (Bard et al., ll999tlBond etail 12001 
iKerrl I200U ISteinhilber et all 120091) Figure 6 shows that the 
J/S conjunctions in AD 1 146 and 1206 occurred approximately 
at the same position of the conjunction in 2000. Thus, the 



J/S co njunction line presents a 800-860 year cycle, as Kepler 
(1606) found where the tropical orbital periods (Pj - 11.857 
and Ps = 29.424) were used. Moreover, the distance of the 
two planets from the Sun at the conjunction dates presents a 
quasi 900-year cycle because of their orbital precession, as 
it is evident by comparing the tidal elongations in 1086 and 
2000, which could be obtained using the sidereal orbital peri- 
ods (Pj - 11.862 and Ps = 29.457), which are more phys- 
ical relative to the Sun. The orbital pattern reveals the exis- 



tence in the solar system of quasi-millennial cycles. Quasi- 
millennial cy cles have been found in sol a r proxy reco r ds r fo r 
example se e: Hood and Jirikowic ( 1991 ): Charvatova (2000); 
Bond et alJ d200ll) : iKerj d200lh : ISteinhilber et al.1 d2009l) l and 
quasi 900- year cycles have been fo und in the Holocene climate 
variability ( Schulz and PaulU2002h . A quasi-millennial cycle is 
quite evident in numerous global temperature reconstructions 
that show a large preindustrial variability. These reconstruc- 
tions show a distinctive Roman Optimum (0-400 AD), a Me- 
dieval Dark Age (400-900 AD), Medieval Warm Period (MWP) 
(900-1300 AD), a Little Ice Age (LIA) (1 400-1800 AD) and a 
Current Warm Period (CWP) si n ce 1800 dMoberg et al] 1 2005 



Loehle and Mc Culloch, 2008; McShane and Wyner, 2011 



Liungqvisl " 2010l) and in numerous recent regional records 



such as from England, Sargasso sea, Japan, Antarctic, Canada, 
Chile, Indo-pacific region, China, tropical South America, and 

I 1 r-^ 1 I ' *— i I 1 I 1 I 1 

many more (Lamb, 1965; Parker etal., 1992; Keigwin, 1996; 
Adhikari and Kumon, 2001; Khim et al., 2002; Huang et al 
20081 : ICook et all 120091; Ivon Gunten et all [20091: lOppo et al 



2009HGe et alll2010llKellerhals et alll2010l:lMann et alll2008 
Liungqvistl 120091) . However, these quasi-millennial cycles, as 
well as other secular and multi-secular cycles may be the topic 
of another work; we only note that the latter temperature re- 
constructions contradict the Hockey Stick temperature graphs 
(Man n and JonesL 120031) that do not show a significant MWP 
nor a LIA, but only an unprecedented CWP. The most recent 
temperature reconstructions would imply the overall solar forc- 
ing of the climate has been underestimated b y at least a facto r 
of three by traditional energy balance models IScafettal (1201 (3). 
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7. Cloud cover oscillation as a possible mechanism for cli- 
mate cycles 

Total solar irradiance variations alone are usually cl aimed 
to be too small to induce large climate variations (IPCC, 
2007b . Additional solar related forcings would be necessary 
to explain climatic variations. Indeed, there may be an ad- 
ditional effect on climate due to ultraviolet modulation of the 



strato s pheric ozone and strat ospheric water vapor (Stub er et al. 
200 1| ISolomon et all [2010), and a modulation of the cloud 
cover due to cosmic ra y flux variation an d a variation of the 
global electric circuit (Shaviv and Veizerl Eool ISvensmarkL 
l2007tlRohs et alll2010l:lLaken et all 1201a iTinslevl l2008h . 

Variations of the frequency of mid-latitude aurora reasonably 
reveal electric variations occurring in the terrestrial magneto- 
sphere and ionosphere. Likely, when the upper atmosphere is 
more ionized, it is more sensitive to solar wind peak activity and 
mid-latitude aurora events occur more frequently. Therefore, 
the mid-latitude aurora records can be considered as a proxy of 
the variation of the atmospheric ionization and, likely, of the 
atmos pheric g l obal electric circuit. As proposed by several au- 
thors (ITinslevl |2008e ISvensmaria 120071) . changes in the global 
electric circuit can drive changes of cloud formation and cover. 
A change of cloud cover alters the albedo and can potentially 
cause a significant climate change. Because auroras oscillated 
with planetary frequencies the climate would oscillate with the 
same frequencies and would be almost perfectly synchronized 
to the 10-11, 20-21 and 60-62 year cy cles of th e J/S co njunc- 
tion, tidal and solar cycles, as shown in lScafetta d2010bh . Sev- 
eral other cycles are surely present, but we do not discuss them 
here. 

To emphasize the effect that a small variation of the albedo 
may have on the climate sensitivity to solar changes we can do 
the following exercise. Let us write: 



AT 



F 



5> 



AF, 



(4) 



where A7> is the change in temperature caused by a small 
change of a generic forcing AF. The value kp is the equilib- 
rium climate sensitivity to a small changes in the forcing, as 
measured by AF. In the case of solar irradiance the above equa- 
tion is substituted with AT$ ~ k$ AFs ■ The same can be done 
for any other forcings. 

By differentiating directly a corrected Stefan-Boltzmann's 
black-body equation 



we get 



AT T , 
k s = — = — = 0.053 K/Wm- 2 
AI 41 



(5) 



(6) 



using / = 1360 W/m 2 as the average total solar irradiance, 
a = 0.3 as the average albedo, T - 289 K as the global average 
temperature of the Earth's surface and cr = 5.67 x 10~ 8 W/m 2 K 4 
is the Stefan-Boltzmann's constant. A corrective amplification 



coefficient / = 1 .65 is used to approximate the greenhouse ef- 
fect of the terrestrial atmosphere to let Eq. |5]to give an equilib- 
rium temperature of T — 289 K. 

The above ks estimate in Eq. [6] assumes, for example, that 
the albedo a is constant. Now, let us assume that an increase of 
the total solar irradiance input of AI = 1 W/m 2 is accompanied 
at equilibrium with a decrease of the albedo by just 1% that 
could be caused also by alternative related solar forcing such as 
a cosmic ray flux solar modulation. The decrease of the albedo 
may be due to a decrease of the low cloud cover and other pos- 
sible factors such as a darkening of the Earth's surface due to 
a melting of the glaciers and biosphere changes. If we assume 
that everything else remains constant, the climate sensitivity to 
total solar irradiance change becomes 



k = — ~ 

S AI 



4/ (l-fl«0.99)(/+A/)/ _ if 

~AI~ 



(l-a)// 
4cr 



0.36 K/Wm- 2 , (7) 



which is seven times larger than the value found in Eq. [6] and 
would be approximately compatible with the empirical values 
found in Scafetta (2009). This suggests that even very small 
solar-induced modulation of the albedo (just a 1 or 2% varia- 
tion) can greatly amplify the climate sensitivity to solar irradi- 
ance changes even by one order of magnitude and could eas- 
ily induce climatic oscillation of the order of a few tenths of 
Celsius degree, which would be perfectly compatible with the 
observations. 

There may also be the possibility that the terrestrial albedo is 
directly modulated by other astronomical mechanisms such as 
a compression of the Earth's magnetosphere due to tidal forces. 
Also in this case Eq. (0 would present an albedo that would 
be a function of time, a(t), and this function would vary with 
the periods of the astronomical forcing. Even in the eventuality 
that the incoming total solar irradiance / and the greenhouse 
feedback coefficient / remain constant, it is easy to prove, using 
Eq. ©, that a reduction of the albedo by just 1 %, for example a 
decrease from a = 0.300 to a — 0.297, would cause an increase 
of the global temperature at equilibrium of about 0.3 K. 

The oscillations observed in the aurora records may suggest 
that the cloud cover oscillates in the same way. Thus, we can 
assume, in first approximation, an albedo function that changes 
periodically as 



a(t) — fl() + a i cos ( w i f + 0i) i 



(8) 



where oq ~ 0.3 is the average terrestrial albedo, <y, and 0, with 
i = 1, 2, 3, ... are the planetary frequencies and phases, and a,- 
are the amplitudes of these cyclical fluctuations. 

We observe that a significant, although small decadal modu- 
lation, of the low cloud cover since 1980, appears synchronized 
to the decadal solar cyc l e dHarrison and Stephenson! |2005; 
Svensmarkl l2007t iBrowni l2008h. Also a ca reful analysis of 
a critical study ( Sloan and Wolfendalel 2008) would approxi- 
mately support the theory: their figure 1 shows the low cloud 
cover anomalies in the polar, intermediate and equatorial re- 
gions from 1980 to 2005, and a sufficiently clear decadal cycle, 
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Figure 9: Connections between weather and climate with thunderstorms, solar 
activity, and galactic cosmic ray flux, via the global atmospheric electric circuit 
and cloud and aerosol microphysics. Solar and heliospheric activity is supposed 
to be partially driven by the gravitational and magnetic fields of the planets, as 
they move around the Sun. See Tinsley (2008) for more details about the single 
mechanisms reported in the figure. 



in synchrony with the decadal solar cycle, is observed. Sloan 
and Wolfendale's figure 1 suggests that the decadal cycle in 
cloud cover has amplitude of about 2%, and it is in synchrony 
with the decadal solar cycle. Moreover, the same figure shows 
a downward multidecadal trend in cloud cover of about 4% 
from 1980 to 2000 that could well agree with the 60-year mod- 
ulation of the auroras and in the global temperature observed 
above in Figure 2. The percentages relative to the decadal and 
60-year modulation would approximately correspond to about 
0.5% and 1 % change of the Earth's albedo respectively because 
these clouds may cover about 30-40% of the Earth. Indeed, 
IWildl d2009h showed that a global solar brightening occurred 
from 1910 to 1940 and from 1970 to 2000, and global solar 
dimming occurred from 1940 to 1970: see his figure 9. This 
60-year reversal from solar brightening to dimming and from 
solar dimming to brightening was likely due to a cloud cover 
60-year modulation and would support the theory herein pre- 
sented. 

An about 0.5% and 1% albedo change would greatly amplify 
the climate sensitivity to total solar irradiance, as shown in the 
above exercise and could also be quantitatively compatible with 
a 0.1 °C climate signature associated to the decadal astronomi- 



cal cycle ( Scafetta , 2009) and with a 0.3/0.4 °C natural warming 
observed from 1970 to 2000, as deduced from the 60-year cycli- 
cal modulation of the temperature depicted in Figure 2. These 
results would be compatible wit h the empirical climate sensitiv- 
ity estimated by Scafetta (2009). Thus, the observations appear 
to be consistent with the theory of an oscillating albedo. 

Figure 9 summarizes a chain of mechanisms that, accord- 
ing to the finding of this paper, reasonably links the plane- 
tary motion to climate change through solar and heliosphere 
modulation of the magnetosphere and ionosphere that regulate 
the cloud cove r percentage. The figure updates a figure from 
Tinslevl Q2008) where more details about the single physical 
mechanisms are found. The planets, in particular Jupiter and 
Saturn, may modulate solar activity, the heliosphere and the 
magnetosphere in proximity of the Earth. Therefore, a cyclical 
modulation of the electric properties of the ionosphere occurs, 
as revealed by the aurora record. This modulation causes a syn- 
chronized change of the cloud system through the atmospheric 
electric circuit. The modulation of the cloud system, together 
with the solar variations, drives all other climate subsystems 
inclu ding the ocean. All subsystem s of the climate synchro- 
nize (Scafetta, 2010b; Strogatz, 2009). The resulting output is a 
global climate system that oscillates with the planetary frequen- 
cies. An additional modulation with a period of about 9.1 years 
is caused by the solar/lunar tidal forces that would act more di- 
rectly in the ocean and in minor degree on the magnetosphere 
and ionosphere. 

In addition, the variation of the length of the day 
(LOD) of the Earth has been proposed as another mech- 
anism linking planetary m otion to climate change. LOD 



presents a 6 0-yea r cycle ([S tephenson and Morrison 
KlvashtorinL 120011 IRoberts et all 120071: IMazzarella . 



1995 



2008 



Mazzarel iaand Scafettall201 l[lMorner[l2010HScafettal[2010bh . 
Because the LOD 60-year modulation is negatively correlated 
to the 60-year climate cycle the following mechanism may be 
working. When the low cloud cover decreases there is an in- 
crease of the energy stored in the oceans that causes variations 
in atmospheric zonal wind circulation and eventually causes 
more water to move from the equator to the poles, causing 
a global warming. This toward-pole mass movement induces 
a LOD decrease because the Earth, by becoming a little bit 
more spherical, speeds up because of the decrease of its rota- 
tional inertia. The latitudinal distribution and transport of en- 
ergy and momentum associated to changes of zonal winds in- 
duced by cloud cover solar-induced changes have been recently 
advocated also for explaining the presence of a significant so- 
lar sig nature in the semiannual LOD variation ( Le Mouel et all 
201(£ 



8. Reconstruction and forecast of the climate oscillations 
based on aurora cycles 

According Eq. [4] if we can assume that solar irradiance 
and/or the albedo periodically oscillate, then the temperature 
at the surface would also oscillate with a similar set of frequen- 
cies. In fact, in first approximation, the surface temperature 
variation would be proportional to the incoming irradiance AF. 
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However, the same result would be obtained also if the heat ca- 
pacity of the Earth is taken into account. For example, a simple 
energy balance model could be given by a first order differential 
equation of the type: 



dAT(t) _ kAF(t) - AT(t) 
dt t 



(9) 



where A; is a climate sensitivity to th e forcing and r is the ther- 
mal relaxation time of the system (IScafettaL 120091) . As well 
known, if the forcing is of the type AF(f) = cos(wf), the solu- 
tion of the above equation is 



AT(t) = k 



cos(wf) + wt sin(<yf) 

1 + L0 2 T 2 



+ c - A cos( cot + cb) + c(10) 



where ct and A are given phase and amplitude that depend on 
co and r, and c is a given constant. Indeed, this result emerges 
from quite general mathematical properties under the approxi- 
mation of small forced oscillations, and would be also obtained 
with professional climate GCMs. Thus, if the solar irradiance 
and/or the albedo are oscillating with a given set of frequen- 
cies, then the surface temperature too would present a harmonic 
component made of the same set of frequencies (in general 
there may be also sub- and super-harmonics). We would have: 



AT{t) = A + V A; cos(w,f + <pi). 



(ID 



Therefor e, as already note d by Lord Kelvin to solve the tidal 
problem ( Thomson! 1881), if the purpose is only to determine 
the response of the climate system to a harmonic forcing, gen- 
eral mathematical theorems allow us to bypass the problem of 
accurately solving the single physical mechanisms. Indeed, the 
parameters A, and <pi, which would contain all physical infor- 
mation about the climate system, could be easily measured and 
determined by simple regression on the actual temperature ob- 
servations. 

It is necessary to test whether a harmonic model based on 
aurora cycles can efficiently both reconstruct and forecast the 
climate oscillations. In fact, a regression model that uses a set of 
parameters would be able to fit only the region used to calibrate 
the model, but if the model is not made of the right dynamical 
components it will fail any forecasting. 

We proceed by assuming, in first approximation, that the 
global surface temperature record has a harmonic component 
characterized by a basic set of five periodic signals with periods 
equal to 9.1, 10.5, 20, 30 and 60 ye ars, as found in the aurora 
records and in the planetary motion (IScafet ta. 2010b). Because 
the temperature rose from 1850 to 2010, this upward secular 
trend needs to be detrended from the data because the proposed 
harmonic model is limited to decadal and multidecadal cycles. 
We use two fit functions of the type: 



Fi(f) = Ai(r- 1850) 2 + fli 



and 



F 2 (?) = A 2 (f- 1850) 2 + fi 2 



(12) 



(13) 



We use Fi(t) to fit the temperature record from 1850 to 2010. 
We use F2(t) to fit the temperature record from 1850 to 1950. 



We find that A x = 2.9 ■ 10~ 5 + 2 ■ 10~ 6 "C/y 2 , B x = 
-0.41 ± 0.02 °C, and that A 2 = 2.7 • 10~ 5 ± 2 ■ 10~ 6 °C/y 2 , 
B 2 = -0.40 + 0.02 °C. We observe that the two acceleration 
coefficients a\ and a 2 do not drastically change. This fact may 
indicate that in 1950 a quadratic fit of the temperature from 
1850 to 1950 would have well forecasted the observed 1950- 
2010 temperature upward trend, but this result is likely a coin- 
cidence. 

Indeed, the secular upward warming trend may be 
due to a combination of a millenarian cycle, of a bi- 
secular cycle, of poorly corrected urban heat is land effects 
dMcKitrick and Michael! 120071: lMcKitric]d.l201Ch and o f an- 
thropoge nic GHG emissi ons. IScafetta! d2009i l2010al) and 
Loehle and Scafettal (1201 lb showed that at least 50% of the up- 
ward secular trend of the temperature since 1900 could be in- 
duced by the correspondent secular increase of the solar activity 
during this same period. Here we do not investigate further this 
issue because the above quadratic polynomial fit is not part of 
the astronomical model, but only a convenient way to represent 
the trending of the temperature from 1850 to 2010, not outside 
this period. 

That the harmonic model may likely work is already implicit 
in Figure 10. Here, it is observed a very good synchrony be- 
tween the quasi 20 and 60-year oscillations of the solar system, 
which are mostly driven by Jupiter and Saturn, and the corre- 
spondent oscillations observed in the climate system. 

The astronomical harmonic model must be tested on its fore- 
casting capabilities to be credible. We consider two indepen- 
dent time intervals (1850-1950 and 1950-2010): 



lln(t-T{)\ lint \ 

Pi(t) = oicosl — I + bi cos I — + ari I + 

lln{t-T x )\ lint \ 

Ci COS + d\ COS hOl + 

\ 20 ) \IQ.5 H ) 

tint \ 
eicosl — +yi|+/i 



and 



lln(t-T 2 )\ , tint \ 
P 2 (t) = a 2 cosl — l + /? 2 cosl — + a 2 \ + 

lln(t-T 2 )\ lint \ 

c 2 cos + d 2 cos 1- B 2 + 

\ 20 ) \10.5 H ) 

tint \ 
<?2COS I — +y 2 \+f 2 



In our simplified model, we use T\ = 2000.5 and T 2 = 1941, 
which correspond to the dates of the last two J/S conjunctions 
of Group I, see Table 2. This choice is also supported by the 
finding depicted in Figure 10. 

We proceed by analyzing three periods. We use the func- 
tion F\ (t) + P\(t) to fit the temperature data during two periods: 
from 1850 to 2010 and from 1950 to 2010. We use the func- 
tion F 2 (t) + P 2 (t) to fit the temperature data from 1850 to 1950. 
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Figure 10: [A] Rescaled 60-year modulation of the solar speed relative to the 
solar system barycenter (SSSB) (black) (see Figure 6B) against the global sur- 
face temperature record detrended of its quadratic fit. [B] Rescaled modulation 
of the solar speed relative to the solar system barycenter (black) (see Figure 6B) 
against two alternative pass-band filtered records of the temperature around its 
two decadal oscillations. The figures clearly indicate a strong coherence be- 
tween the astronomical oscillations and the oscillations observed in the climate 
system. 



The coefficients of the three regression models are reported in 
Table 3. We find that the three sets of regression coefficients 
coincide within the error of measure: both the amplitudes and 
the phases coincide. The reduced^ is 0.34 with 9 degrees of 
freedom and the coherence probability between the two mod- 
els is Pg(x 2 > Xo) > 95%. Thus, within the limits of the data 
herein used, it is found that the astronomical harmonic model 
can both reconstruct and forecast the temperature signal within 
a high statistical confidence level. 

Figure 11A depicts the three regression model outputs 
against the global surface temperature, which has been 
smoothed with a four-year moving average algorithm. Figure 
11B depicts only the oscillating functions Pi(t), which fits the 
periods 1850-2010 (solid) and 1950-2010 (dash), and the func- 
tion P 2 (i), which fits the period 1850-1950 (dots). Figure 11A 
and 11B clearly show a very good matching during the entire 
period 1850-2010 between the decadal and multidecadal cli- 
mate oscillations and the oscillations reconstructed by the mod- 
els in all three cases. 



Figure 1 1 : Astronomical harmonic constituent model reconstruction and fore- 
cast of the global surface temperature. [A] Four year moving average of the 
global surface temperature against the climate reconstructions obtained by us- 
ing the function F\(t) + Pi(r) to fit the period 1850-2010 (black solid) and 
the period 1950-2010 (dash), and the function F 2 {t) + P%(t) to fit the period 
1850-1950 (dots). [B] The functions Pi{t) and P 2 (.t) represent the periodic 
modulation of the temperature reproduced by the celestial model based on the 
five aurora major decadal and multidecadal frequencies. The arrows indicate 
the local decadal maxima where the good matching between the data patterns 
and the models is observed. Note that in both figures the three model curves 
almost coincide for more than 200 years and well reconstruct and forecast the 
temperature oscillations. 



The result depicted in the figure implies that if we were in 
1950, by using an astronomical harmonic constituent model 
based on aurora cycles and on the temperature record from 1850 
to 1950, it could have been possible to forecast with a reason- 
able accuracy the climate oscillations within a decadal and mul- 
tidecadal scale occurring from 1950 to 2010. The figure also 
implies the opposite, that is, that the decadal and multidecadal 
climate oscillations from 1850 to 1950 could be forecasted by 
using the same astronomical harmonic constituent model cali- 
brated on the temperature data from 1950 to 2010. Note that 
during the period 1915-1940 the decadal cycle appear to dis- 
appear because during this period the 9.1 and 10.5 year cycles 
interfere destructively. The remarkable result depicted in Fig- 
ure 1 1 does imply that the climate oscillations within the secular 
scale can be accurately forecasted by using the major decadal 
frequencies found in the mid-latitude aurora records, which ap- 
pear to be driven by astronomical cycles associated to the move- 
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ment of Jupiter and Saturn (about 10, 20, 30 and 60 years), to 
the lunar cycle (9.1 years) and to the 1 1-year solar cycle. 

About the imminent future, Figure 11A suggests that if the 
quadratic background upward trend of the temperature con- 
tinues during the next decades, the global temperature may 
remain approximately constant until 2030-2040, and this re- 
sult would be approximat ely compatible with that found in 
Loehle and Scafettal ( 201 1 ) by using a different methodology 



that yield a net anthropogenic warming component of about 
+0.66 °C/century from 1940 to 2050 against the IPCC esti- 
mate of +2.3 "CI century during the same period. However, the 
quadratic trend used here should not be trusted after 2040 also 
because is not part of the harmonic model, but only a convenient 
way to fit the temperature from 1850 to 2010. 

The observed warming since 1850 could have been partially 
caused by a large quasi-millennial climate cycle, which has 
caused the Medieval Warm Period, the Little Ice Age and the 
Modem Warm Period. Probably this large quasi-millennial 
cycle just entered or is going to enter into its cooling phase 
and, therefore, it can further cool the planet on a multi-secular 
scale. Other cooling natural component can derive from a bi- 
secular natural solar induced cycle. These cycles are known 
to have regulated multidecadal cold periods (the well-known 
P ort, Wolf , Sporer , Maunder, Dalton minima). In fact, as noted 



Scafetta (2010a) the last four sunspot cycles (#20-21-22-23) 



look quite similar to the four sunspot cycles (#1-2-3-4) that oc- 
curred just before the Dalton minimum. That the Sun may be 
enter ing in a prolonged period o f minimum activity is reason- 
able (iDuhau and de Jagei . 2010h . That a multi-secular climate 
natural cycle may be turning down is also indirectly suggested 
by the fact that since 1930 th e sea-level accelera t ion ha s been, 
on average, slightly negative ([Houston and Deanl 1201 lb . Inter- 
estingly, the measured deceleration of the sea-level rise from 
1930 to 2010 occurred despite the strong positive accelera- 
tion of anthropogenic GHG emissions during the same period, 
which would indicate that other more powerful forcings (such 
as an astronomical forcing of the cloud system) are the major 
regulators of climate change. 



9. Conclusion 

Four centuries ago, Johannes Kepler explained that earthly 
nature couldn't help but respond to the dictates of heavenly 
harmonies, and said that nature is affecte d by a n aspe ct just 
as a farmer is moved by music to dance (IKempL 120091) . Ke- 
pler clearly shared the common belief of his time that the cli- 
mate was influenced by astronomical cycles and understood 
the subtle phenomenon of collective synchronization that has 
been extensively studie d in nonlinear complex science since 
the times of Hu ygens dPikovskv et al. , 2001 : Strogatz , 2009; 



I Scafettal l2010bl) . Indeed, climate change records present ge 
ometrical characteristics that suggest that climate is synchro- 
nized, probably through the Sun and the heliosphere, to com- 
plex astronomical cycles driven by planetary harmonics. 

In this paper, we have studied the historical record of mid- 
latitude auroras from 1700 to 1900, and of the Faroes Islands 



from 1873 to 1966 to search for a possible physical mech- 
anism linking planetary motion to climate. We have shown 
that mid-latitude aurora records and the global surface temper- 
ature record share a set of oscillations with periods of about 
9.1, 10-11, 20-21, 30 and 60-62 years. Other shorter and 
longer oscillations may be present, but they are not discussed 
here. In particular, clear quasi 60-year cycles in the aurora 
records are synchronized to the 60-year cycle observed in the 
global surface temperature and in multi-secular proxy climate 
reconstructions of both the Atlantic and Pacific climatic os- 
cillations since 1650 and in Indian summer monsoon records. 
Charvatova-Jakubcova et al. dl988h found that a large 60-year 



cycle is present in the mid-latitude aurora record also for the 
longer period from 1001 to 1900 and other planetary frequen- 
cies are present as well in the millennial aurora record. Nu- 
merous other studies have found quasi 10, 20 and 60-year cy- 
cles in multiple climate records and in astronomical records, as 
summari zed in the Introd uction. By taking into account the re- 
sults of lScafettal d2010bb . this synchrony exists also with the 
global ocean and land global surface temperature records of 
both hemispheres. 

The aurora record cycles reveal a direct or indirect plane- 
tary influence on the Sun and on the Earths magnetosphere and 
ionosphere. Indeed, proxy reconstructions may suggest a 60- 
year cycle in the total solar irradiance (TSI) was almost stable 
or slightly decreased from 1880 to 1910, and it increased from 
1910 to 1940. From 1940 to 1970 TSI may have de creased as 
the TSI reconstruction of Hovt and Schattenl ( 1997 ) suggests. 
Hoyt and Schatten's TSI reconstruction well c orrelates with 
the temperature records during the last century (ISoonl 120091) . 
Finally, an increase of the solar activity from 1970 to 2000 
and a decrease afterward would be supported by the ACRIM 
TSI satellite composite, whi ch may more faithfully repro - 
duce the satellite obs e rvatio ns Willson and Mordvinov (2003); 
Scafetta and Willsonl J2009I); IScafettal d201 lb (but not by the 



PMOD composite Frohlich and LeanJ dl998l) . which is the TSI 
record preferred by the I PCCl(l2007l) ). Thus, solar activity could 
have been characterized by a quasi 60-year modulation super- 
impo sed to other larger secular, bi-s e cular and millennial cy- 
cles (Kf ivskvl 1 1984c lOgurtsov et all 120021) . although it may 
not appear evident in every total solar ir radiance reconstruction 
dKrivova et al.U2007l:IWang et alll2005l) . 

It is possible that when Jupiter and Saturn are closer 
to the Sun, there may be an increased solar activity be- 
cause of the stronger pl anetary tides and other mechanisms 
dWolff and Patronell2010l) . and a stronger magnetic field within 
the inner region of the solar system form, although the patterns 
may be more complicated because of the presence of other cy- 
cles that will be discussed in another paper. A stronger solar or 
heliospheric magnetic field better screens galactic cosmic ray 
fluxes. Fewer cosmic rays reaching the Earth imply a weaker 
ionization of the upper atmosphere. As a side-effect, less au- 
roras form in the middle latitudes because a stronger magnetic 
field and a less ionized ionosphere mostly constrains the auro- 
ras in the polar region. In addition, the level of ionization of 
the atmosphere has been proposed as an import ant mechanism 
that can modulate the low cloud cover formation ( Tinslevl 2008 
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Kirkbv, 2007; Svensmark et al. , 2009). Essentially, when the 
ionization is weaker, less clouds form. A solar and heliospheric 
modulation of the cloud system would greatly contribute to cli- 
mate change through an albedo modulation (see Eq. [8}. The 
above sequence of mechanisms would explain why the climate 
presents oscillations at multiple frequencies that are synchro- 
nized with the aurora and the planetary cycles. 

We have used a phenomenological harmonic model based on 
five decadal and multidecadal frequencies with periods of 9.1, 
10.5, 20, 30 and 60 years that has been detected in the aurora 
records and that could be associated to evident astronomical 
and solar/lunar tidal cycles. We have shown that it is possible to 
forecast the climate oscillations occurred from 1950 to 201 us- 
ing the climatic information derived from the period 1 850- 1 950 
and the frequency information deduced from the mid-latitude 
aurora before 1900. Analogously, we have shown that it is pos- 
sible to forecast the climate oscillations that occurred backward 
from 1850 to 1950 using the information derived from the pe- 
riod 1950-2010. Thus, these findings strongly support the the- 
sis that the climate oscillations can be approximately forecasted 
by using astronomical cycles. The proposed astronomical har- 
monic constituent model for climate change based on aurora 
cycles is conceptually equivalent to the commonly used tide- 
predicting machines based on planetary harmonic constituent 
analysis conceived by Lord Kelvin in 1867. 

Interestingly, the traditional Chinese, Tamil and Tibe t an cal - 
endars are arranged in major 60-year cycles (Aslaksen, 1999). 
Perhaps, these sexagenarian cyclical calendars were inspired 
by climatic and astronomical observations and were used for 
timing and regulating human business. For example, even an- 
cient SJanskrittexts report about a 60-year monsoon rainfall cy- 
cle Jlvengari r2009) and associate it to the movement of Jupiter 
and Saturn, the Brihaspati 60-year cycle, which may explain 
why Asian populations used sexagesimal calendars. Indeed, a 
60-year cycle linked to Jupiter and Satu rn was e x treme ly well 
known to several ancient civilizations (Temple, 1998). The 
major cycles discussed in this paper are also approximately 
found in the major busines s cycles ( Pustilnik and Din[ 2004; 
Korotavev and Tsirel 2010h . A link between planetary motion, 
climate and economy (which mostly in the past could be driven 
by agricultural productivity) would ultimately explain the in- 
terest of the ancient civilizations in tracking the position of the 
planets and their attempts in developing multiscale cyc l ical ca l- 
endars ( Ptolemy , 2nd century ; Ma 'sari 886; Swerd low , 1998 ). 

In conclusion, the resul ts presented here strongly support and 
reinforce the argument of lScafetta (2010a b) that the climate is 
forced by astronomical oscillations related to the Sun, the Moon 
and the planets, and, as Figure 1 1 shows, a significant compo- 
nent of it can be forecasted within an acceptable uncertainty 
with appropriate harmonic models. 
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Cycle 


Temp/Fig4A 


AAR 


fig4B 


fig4C 


fig4D 


#1 


9.1+0.2 


9.15+0.2 


9.1+0.2 


9.2 + 0.2 


n/a 


#2 


10.4 ±0.3 


10.33 + 0.3 


10.0 + 0.3 


10.2 + 0.3 


10.8+0.3 


#3 


20.9 ± 0.9 


20.6 ± 0.9 


21.3+0.9 


20.3 + 0.9 


20.2 + 0.9 


#4 


32 + 2 


29.5 + 2 


29 + 2 


30 + 2 


n/a 


#5 


62 + 5 


62 + 5 


62 + 5 


64 + 8 


60 + 5 



Table 1: First column: periods of the cycles #1, #2, #3, #4 and #5 found in the global surface temperature: see Figure 4A. Third column: periods of the cycles #1, 
#2, #3, #4 and #5 found in the 1700-1966 aurora record: see Figure 4B. Fourth column: periods of the cycles #1, #2, #3, #4 and #5 found in the 1700-1880 aurora 
record: see Figure 4C. Fifth column: periods of the cycles #2, #3 and #5 found in the 1872-1966 aurora record: see Figure 4C. Second column: average of the 
periods of the cycles #1, #2, #3, #4 and #5 found in the three aurora records. The results depicted in the first and second column are perfectly coherent within the 
error of measure for each couple of values. For the 5 frequency couples combined, the reduced^- 2 is xl = 0.18 with 5 degrees of freedom (P 5 (x 2 > xl) > 96%). 



Group I 


Jupiter 
Right asc. 


Decl. 


Dist 


Saturn 
Right asc. 


Decl. 


Dist 


tidal 
elong. 


Jun 23, 2000 


3h 19m 29.7s 


17° 19' 48" 


4.996 


3h 20m 46.9s 


16" 09' 35" 
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Nov 15, 1940 
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17° 39' 14" 
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Jul 26, 1 146 
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Group II 
















Apr 17, 1981 
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5.451 
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-0°38' 21" 
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Aug 23, 1921 
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May 8, 1802 
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9h 57m 07.1s 


13° 33' 53" 


5.378 
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13" 58' 48" 
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Feb 2, 1683 


9h 18m 23.2s 


16° 36' 37" 


5.347 


9h 18m 43.6s 


16° 50' 56" 


9.177 
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Group III 
















Apr 16, 1960 


19h42m 18.9s 


-21° 40' 29" 


5.146 


19h 42m 10.4s 


-21° 23' 03" 


10.028 


1635 


Sep 28, 1901 


19h 07m 48.3s 


-22° 40' 42" 
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19h 07m 38.6s 


-22° 13' 36" 


10.053 


1607 


Mar 12, 1842 


18h33m 28.3s 


-23° 12' 38" 


5.207 


18h33m 23.7s 


-22" 36' 25" 


10.072 


1580 


Aug 26, 1782 


17h 59m 51.0s 


-23° 16' 44" 


5.236 


17h 59m 50.7s 


-22° 32' 18" 


10.082 


1555 


Feb 11, 1723 


17h26m 46.3s 


-22° 54' 28" 


5.264 


17h 27m 03.5s 


-22° 02' 54" 


10.086 


1531 


Jul 31, 1663 


16h 54m 47.2s 


-22° 08' 09" 


5.291 


16h 55m 15.1s 


-21° 10' 25" 


10.083 


1509 



Table 2: Approximates dates of the J/S conjunctions relative to the Sun, which slightly differ from the conjunction dates as relative to the Earth. The coordinate of 
the planets are in Right Ascension and Declination in heliocentric equatorial coordinates which are relative to the equatorial plane of the Earth. The last column 
report the approximate value of the tidal elongation in Km (see Eq. |2j at the orbit of the Earth, that is, at 1 AU from the Sun. 





Pi (f)(1850-2009) 


Pi(t) (1950-2009) 


P 2 (t) (1850-1950) 


a 


0.103 + 0.011 


0.109 + 0.018 


0.098 + 0.015 


b 


0.022 + 0.011 


0.031+0.018 


0.020 + 0.015 


c 


0.040 + 0.011 


0.043 + 0.018 


0.036 + 0.015 


d 


0.027 + 0.011 


0.030 + 0.018 


0.027 ± 0.015 


e 


-0.053 + 0.011 


-0.055 + 0.018 


-0.050 ± 0.015 


f 


-0.006 + 0.008 


-0.024 + 0.013 


-0.006 ± 0.010 


a 


2.08 + 0.5 


1.84 + 0.6 


2.41 +0.7 




1.43 + 0.4 


1.81+0.6 


1.14 + 0.6 


7 


-0.32 + 0.2 


-0.47 ± 0.3 


-0.26 ± 0.3 



Table 3: Multi nonlinear regression coefficients to fit of the global surface temperature record using the functions Pi(f) and P2(t). Pl(f) is used to fit two periods: 
1850-2009 and 1950-2009. Pi(t) is used to fit the period 1850-1950. Note that the three sets of coefficients coincide within the errors of measure. Thus, there exists 
an accurate statistical coherence between the three models. For the 9 coefficient couples relative to the independent results reported in the last two columns the 
reduced^ 2 is^ 2 = 0.34 with 9 degrees of freedom (Pg(x 2 > xl) > 95%). 
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